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Abstract 

The magnetic properties of Na2CuP207 were investigated by means of ^^P nuclear magnetic 
resonance (NMR), magnetic susceptibility, and heat capacity measurements. We report the ^^P 
NMR shift, the spin-lattice (1/Ti), and spin-spin (I/T2) relaxation-rate data as a function of 
temperature T. The temperature dependence of the NMR shift K{T) is well described by the 
S = 1/2 square lattice Heisenberg antiferromagnetic (HAF) model with an intraplanar exchange 
of J/kB ~ (18 ± 2) K and a hyperfine coupling A = (3533 ± 185) Oe/^lB■ The ^^P NMR spectrum 
was found to broaden abruptly below T ~ 10 K signifying some kind of transition. However, no 
anomaly was noticed in the bulk susceptibility data down to 1.8 K. The heat capacity appears to 
have a weak maximum around 10 K. With decrease in temperatures, the spin-lattice relaxation rate 
1/Ti decreases monotonically and appears to agree well with the high temperature series expansion 
expression for a S =1/2 2D square lattice. 

PACS numbers: 75.10.Pq, 75.40.Cx, 76.60.-k, 76.60. Cq 
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I. INTRODUCTION 



Low-dimensional spin systems with antiferromagnetic interactions have received consid- 
erable attention because of pronounced quantum mechanical effects which result in magnetic 
properties which are quite different from those of the three-dimensional (3D) antiferromag- 
netic substances. According to the Hohenberg-Mermin- Wagner theorenti, one-dimensional 
(ID) and two-dimensional (2D) spin systems with Heisenberg interaction and finite-range 
coupling between spins cannot have long-range order (LRO) at any temperature different 
from zero. In ID systems it does not occur even at zero temperature. However for 2D 
systems at zero temperature, LRO is not forbidden by this theorem. In fact, LRO has been 
rigorously established for spin 5* > 1/2, on a square lattice with nearest-neighbor coupling.i^ 
For the case of 5* = 1/2, there is no solid proof but there are strong theoretical arguments 
that LRO exists."^-^ 

The case of a = 1/2 system on a square lattice, with nearest-neighbor antiferromag- 
netic coupling, has been of special interest because of its proximity to the high-Tc cuprates. 
The interest in quasi- ID spin systems has been stimulated by the hope that better under- 
standing of ID systems might lead to insights into the 2D systems and high-temperature 
superconductivity. 

So far, a large number of = 1/2 compounds have been experimentally investigated, 
which could effectively be described by 2D HAF models. Among the important ones 
are cuprate compounds such as La2Cu04 and YBa2Cu306 which have CuO planes.-^^ In 
La2Cu04, the intraplanar exchange coupling J/ks was reported to be 1800 K. In spite of 
a small interplanar coupling J 10^^ J, the large correlations in CuO planes lead to 3D 
magnetic order at the Neel temperature T/v ~ 300 K. In the above cases T/v is strongly 
dependent on J/ J ratio. 

Most recently, KCuFs has been experimentally investigated as one of the quasi-lD HAF 
compounds.-S^ Unfortunately, this material has a relatively large coupling ratio J/ J ~ 
1.0 X 10-2, due to which, the Tn / J ratio (~ 39 K/203 K) was found to be large. Later 
Sr2Cu03 was reported to be another ID HAF system,— with a significantly reduced T/v/ J 
ratio of ~ 5 K/2200 K ~ 2 x 10^^. Due to the relatively small Tjv value compared to 
the exchange coupling J/ks, ID behavior is observed over a wide range of temperature. 
Recently, Nath et alM have reported that Sr2Cu(P04)2 and Ba2Cu(P04)2 are two ID HAF 
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systems, which do not appear to undergo Neel ordering even at a very low temperature 
(T ~ 0.02 K). Their exchange couphng constants have been reported to be 165 K, and 151 
K respectively. It is thus of interest to synthesize and characterize additional S = 1/2 ID 
or 2D HAF compounds to improve our understanding of such systems. In this paper, we 
present a detailed study of the magnetic properties of Na2CuP207 via susceptibility and 
^^P nuclear magnetic resonance (NMR) experiments. Our objective is to report detailed 
magnetic measurements on a new potentially low-dimensional system and to analyse the 
data based on available models. An additional objective is then to motivate the theorists 
to model real (and more complex) systems as the ones we report and compare the result 
of their simulations with our data. Our results indicate that the magnetic properties of 
Na2CuP2 07 are consistent with the 2D square-lattice S = 1/2 HAF model somewhat better 
than the ID model. In the next section, an overview of the schematic structure is presented 
which motivated us to work on this system. This is followed by details of our experiments. 
The "experimental details" section is followed by results of our magnetic and heat capacity 
measurements accompanied by an analysis and a conclusion. 

II. STRUCTURE 

The structural properties of Na2CuP207 have been reported by Etheredge et al}^ and Er- 
ragh et alM It (the high-temperature phase) crystallizes in a monoclinic unit cell with space 
group C2/c- The reported lattice constants are 14.715 A, 5.704 A, and 8.066 Arespectively 
along a-, b- and c- directions. From the schematic diagram of the structure shown in Fig. 1 
(a) and (b), it is seen that the exchange interaction between Cu^'*" ions could arise due to 
two interaction paths: (i) each CuOg octahedron shares its corners with two similar kinds 
of PO4 groups. The corner sharing takes place in one direction forming [Cu(P04)2]oo chains 
along the c-direction and in this case the magnetic properties would be those of a ID HAF 
chain, (ii) alternatively, there exist (nearly) 180° Cu-O-Cu linkages in the be plane. De- 
pending on the orientation of the d^-^-y-i orbitals {i.e. whether perpendicular to the he plane 
or in the he plane), the system will behave either as a chain-like or a planar magnetic sys- 
tem. In the latter case however, one should note that the Cu^"^ ions are arranged in a face 
centered manner (see Fig. 1(b)). One can, therefore, think of the he plane as comprising 
of two sub-planes (shown in Fig. 1(b) by thick and thin lines). While the Cui— O— Cui 
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(or Cu2— O— CU2) bond angle is almost 180°, the Cui— O— Cu2 bond angle is somewhat less 
than 90°. The deviation of the magnetic properties of Na2CuP207 from the square planar 
case would depend on the relative strength of the Cui— O— Cu2 interaction with respect to 
that of the Cui— O— Cui (or Cu2— O— CU2) interaction. Further, the interaction between 
the be planes is expected to be weak since the distance between them is about 8 Awhich is 
nearly twice the Cu^"*"— Cu^"*" intraplanar distance of 4 A. Also, unlike the 180° Cui— O— Cui 
intraplanar bonds, there appear to be no similar interaction paths perpendicular to the be 
plane. 

III. EXPERIMENTAL DETAILS 

Polycrystallinc sample of Na2CuP207 was prepared by solid state reaction technique using 
NaH2P04.H20 (98% pure) and CuO (99.99% pure) as starting materials. The stoichiometric 
mixtures were fired at 800 °C for 120 hours, in air, with several intermediate grindings and 
pelletization. Formation of nearly single phase sample was confirmed from x-ray diffraction, 
which was performed with a Philips Xpert-Pro powder diffractometer. A Cu target was 
used in the diffractometer with Xav = 1.54182 A. An impurity phase was identified to 
be NasPsOg and the intensity ratio {limp/ 1 max) was found to be 0.05, where limp is the 
intensity of the most intense diffraction peak for NaaPsOg and Imax is that of Na2CuP207. 
Lattice parameters were calculated using a least-square fit procedure. The obtained lattice 
constants are 14.703 (4) A, 5.699 (2) A, and 8.061 (3) A, respectively along a-, b- and c- 
directions. These are in agreement with previously reported values. Magnetization (M) 
data were measured as a function of temperature T (1.8 K < T < 400 K) and apphed 
field H {0 < H < 50 kG) using a SQUID magnetometer (Quantum Design). The heat 
capacity was measured with a PPMS set-up (Quantum Design). The NMR measurements 
were carried out using pulsed NMR techniques on ^^P nuclei (nuclear spin 1 — 1/2 and 
gyromagnetic ratio 7/271 = 17.237 MHz/Tesla) in a temperature range 2 K < T < 300 K 
using a ^He cryostat (Oxford Instruments). We have done the measurements in an applied 
field of about 55 kG, which corresponds to a radio frequency (rf) of about 95 MHz. Spectra 
were obtained by plotting the echo integral (following a 7r/2 — tt pulse sequence with a 7r/2 
pulse of width 4 fis) as a function of the field at a constant frequency of 95 MHz. The 
NMR shift K{T) — [Href — H (T)] /H{T) was determined by measuring the resonance field 



4 



of the sample {H (T)) with respect to a reference H3PO4 solution (resonance field Href)- 
The ^^P nuclear spin-lattice relaxation rate (1/Ti) was determined by the inversion-recovery 
method. Nuclear spin-spin relaxation rate (I/T2) was obtained by measuring the decay of 
the transverse nuclear magnetization with a variable spacing between the tt/2 and the vr 
pulse. 



IV. RESULTS AND DISCUSSION 

We first present the results of our ^^P NMR measurements in Na2CuP207. Since there is 
a unique ^^P site, the ^^P NMR spectra consist of a single spectral line at high temperatures 
(T > 5 K) as is expected for 1=1/2 nuclei. The observed peak position shifts with respect 
to Href in the field-sweep spectra. The temperature dependence of the ^^P NMR shift 
apparently arises due to the temperature dependence of the spin susceptibility XspiniT) via 
a hyperfine coupling to the Cu^"*" ions. The NMR shift is not affected by small amounts of 
extrinsic paramagnetic impurities whereas in the bulk susceptibility they give rise to Curie 
terms. NMR shift data as a function of temperature are shown in Fig. 2, for 5 K < T < 300 
K. They exhibit a broad maximum at 20 K, indicative of short-range ordering. As explained 
earlier, the dominant magnetic behavior of Na2CuP207 could be that of a HAF chain or a 
plane. Consequently, we tested both the 2D (planar) and ID (chain) models to fit the NMR 
shift data. A high temperature (^5^ ^0-7) series expansion for the inverse susceptibility 
l/XspiniT) for 2D S* = 1/2 HAF square lattice was given by Rushbrooke and Wood^^, which 
has the form, 

(1) 



'^/XspiniT) — j^j 2 „2 



4 \n-l 
n=l ^3-^'' 



where x = g is the Lande ^f-factor, /i^ is the Bohr magneton, A^^ is the Avogadro 
number, and C„ are the coefficients listed in Table 1 of Ref. 13. Similarly Johnstoni^ 
parametrized the low-temperature (^^^ < l) simulations of Takahashi^^ and Makivic and 
Dingi^ to obtain 

J 



XspiniT) = ^^^b9 ^0.043669 + 0.039566a; - Q.bUlx^ + 4.684a;^ - 11.13a;^ + lO.SSx*^ - 3.56x^] 



(2) 

For ID HAF chains, the temperature dependence of the susceptibility Xspm (T) was numeri- 
cally calculated by Bonner and Fisher— which accurately predicts the susceptibility for high 
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temperatures (-^ > 0.5). Below, we use the form as given by Estes et alM. 

Ng^nl (0.25 + 0.074975S-1 + 0.075235a;-2) 

^'^'"^ ~ kex ^ (1 + 0.9931X-1 + 0.172135x-2 + 0.757825x-3) ^ ^ 



Since the temperature dependence of the Xspin{T) is reflected in the NMR shift K{T), one 
can determine the exchange coupling J/ks and the hyperfine interaction A simultaneously 
by fitting the temperature dependence of K to the following equation, 

where Ko is the chemical shift. Figs. 2 and 3 show fitting of the ^^P NMR shift data to Eq. 
H taking 

Xspin for HAF square lattice (Eq. ^ and linear chain (Eq. ISJ respectively. In Fig. 
2, the fitting to the 2D high-temperature series expansion was done for 15 K < T < 300 
K whereas in Fig. 3, the experimental data were fitted to the linear chain model in the 
temperature range 5 K < T < 300 K. The parameters extracted from the fit are listed in 
Table I. The Lande g'- value was found to he g = 2.1, which is a typical value for cuprates. 

TABLE I: Values of the parameters obtained from the fitting of NMR shift to Eq. 4 considering 
2D square planar HAF (Eq. 1) and ID HAF (Eq. 3) models. 





Kq (ppm) 


A (Oe/fiB) 




Fig. 2 


100 ± 50 


3533 ± 185 


18 ±2 


Fig. 3 


143 ± 60 


3479 ± 200 


28 ±5 



From Figs. 2 and 3 it is seen that our K{T) data fit somewhat better to the 2D HAF 
model. At low temperature (T < 25 K), the ID fit deviates from the experimental data. 
In Fig. 2, we have also plotted the simulated low temperature curve using Eqs. |21 and E] 
with the parameters Kq, A, and J/ks obtained from the high temperature fit along with our 
experimental data. It is clearly seen that our experimental data do not deviate significantly 
from the simulated curve down to T ~ 10 K while a large deviation is seen below 10 K. This 
suggests some transition or crossover below 10 K. 

An observation of the ^^P NMR lineshapes below 10 K (shown in Fig. 4) reveals a huge 
broadening at lower temperatures. Further, the lineshape develops shoulder-like features 
and finally at 2 K the overall extent of the spectrum is about five times that at 10 K, with 
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at least three distinct peaks. Either a structural or a magnetic transition might be the cause 
for this. No anomaly is seen in the bulk susceptibility (see below) which seems to suggest 
against the occurrence of 3D LRO. 

Magnetic susceptibility x{T) (= M/H) of Na2CuP207 was measured as a function of 
temperature in an applied field of 5 kG (Fig. 5). The amount of ferromagnetic impurity 
present in our sample was estimated from the intercept of M vs. H isotherms at various 
temperatures and was found to be 19 ppm of ferromagnetic Fe^"*" ions. The data in Fig. 5 
have been corrected for these ferromagnetic impurities. As shown in the figure, x{T) exhibit 
a broad maximum at 20 K, indicative of low- dimensional magnetic interactions. With a 
further decrease in temperature, susceptibility increases in a Curie- Weiss manner. This 
possibly comes from defects and extrinsic paramagnetic impurities present in the samples. 
No obvious features associated with LRO are seen for 1.8 K < T < 400 K. 

The broad maximum in x(T) at 20 K could be reproduced by assuming that 

C 

X = XO + + Xspin{T) (5) 

where Xspm{T) is the uniform spin susceptibility for a 5* = 1/2 2D HAF system obtained from 
Eq. H xo is temperature independent and consists of diamagnetism of the core electron shells 
iXcore) and Van-Vleck paramagnetism(xt,^,) of the open shells of the Cu^+ ions present in the 
sample. The Curie- Weiss contribution is ^where C = ^^^^^^^^ ^(hie to paramagnetic 
species in the sample. The parameters were determined by fitting our experimental x{T) 
data to Eq. |3]in the high temperature regime 15 K < T < 400 K. The extracted parameters 
are Xo =(-7 ±2) xlO"^ cmVmole, C =(13±3)xl0-3 cm^K/mole, = 1.7 K, ^ = (18 ±2) 
K, and g = 2.07. Adding the core diamagnetic susceptibility for the individual ions^i (Na^+ 
= -5 X 10"^ cmVmole, Cn^+ = -11 x 10"^ cm^/mole, = -1 x 10"^ cm^mole, O^" = -12 
X 10~^ cm^/mole), the total Xcore was calculated to be —1.07 x 10^^ cm^/mole. The Van- 
Vleck paramagnetic susceptibility estimated by subtracting Xcore from xo is about 3.7 x 10^^ 
cm^/mole, which is comparable to that found in Sr2Cu03 (~ 3.4x10^^ cm^/mole)^^. The 
Curie contributions present in the sample corresponds to a defect spin concentration of 3.5 
% assuming defect spin 5* = 1/2. 

Further, we did heat capacity measurements on Na2CuP207 to look for signs of any 
anomalies at low temperature, signaling a magnetic transition. As seen in the data (Fig. 
6), no sharp peaks are visible which seems to rule out LRO. However, a look at the 



derivative of the specific heat as a function of temperature (see inset of Fig. 6) clearly shows 
a local maximum at about 7 K and a local minimum around 10 K. This suggests that 
the specific heat has an anomaly/peak between 7 K and 10 K.— In the 2D HAF model, a 
broad maximum in the heat capacity is expected at about T = 0.58 J/ (i.e. at about 
10 K in the present case) whereas in the ID HAF model, a broad maximum is expected at 
about T = OASJ/k^ (i.e. at about 16 K in the present case). This further suggests the 
applicability of the 2D HAF model in the present case. 

The time dependence of the longitudinal nuclear magnetization M{t) for ^^P at three 
different temperatures is shown in the inset of Fig. 7. For a spin-1/2 nucleus this recovery 
is expected to follow a single exponential behavior, 

^'i-)-^'it) ^^^J_^),c (6) 



M(oo) "V ^1, 

Our experimental data show good single exponential behaviour over two decades. The spin- 
lattice relaxation rate 1/Ti was extracted from the fitting of the experimental data at various 
temperatures (down to 5 K) to Eq. El- Due to the large line broadening we could not saturate 
the nuclear magnetization below 5 K and hence could not extract reliable 1/Ti below this 
temperature. The temperature dependence of ^^P nuclear spin-lattice relaxation rate thus 
obtained is presented in Fig. 7. With decrease in temperature, it decreases monotonically. 
For 5* = 1/2 2D square lattice, a high temperature series expansion for 1/Ti was given by 
Moriya^^ which has the form, 

Ti^ (1 + J/iAkBT)f'exp [(J/(2fcBT))2(l + J/iAkBT))] 

where, ^ = {^){^^^) taken from Ref. 20. Here, Ah = '^Ah-^ where A is the total 
hyperfine coupling obtained from the experiment. In Moriya's expression a term of the order 
of {J/ksTY occurs in the prefactor and has a negligible effect and hence is not included in 
Eq. |7| Using Eq. [7| and the relevant A and J / ks values obtained from the fit of NMR shift 
to 2D model, we simulated the theoretical curve for Na2CuP207 and it is plotted with our 
experimental data in Fig. 7. Clearly, at high temperatures (T > 25 K), the experimental 
data agree resonably well with the simulated curve. In the case of a 5* = 1/2 ID HAF chain 
model, 1/Ti is expected to be temperature independent at low temperatures (T << J/Zcb) 
due to a dominant contribution from the fluctuations of the staggered susceptibility of the 
ID chain. We are unable to probe this region of temperature due to the large broadening 
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of our NMR spectra seen there. At higher temperatures, from fluctuations of the uniform 
susceptibihty of the ID chain, one expects a hnear variation of 1/71 with temperature. This 
should eventually saturate at even higher temperatures when the spin-susceptibility becomes 
Curie-like. Qualitatively speaking, the observed 1/Ti data could be explained based on the 
above. However, an analytical expression for the temperature dependence of 1/Ti is not 
available in the temperature regime of our experiment and hence no curve-fit is shown in 
the figure. 

The spin-spin relaxation was measured as a function of separation time t between 7r/2 
and TT pulses by monitoring the decay of the transverse magnetization. I/T2 at different 

temperatures was obtained by fitting of the spin-echo decay to the following equation. 



Inset of Fig. 8 shows the spin-echo decays for different temperatures. The extracted spin- 
spin relaxation rates I/T2 arc plotted as a function of temperature in Fig. 7. It is to 
be seen that below about 55 K, the spin-spin relaxation rate I/T2 falls sharply towards 
low temperatures. No indication of 3D LRO was found down to 6 K. The origin of this 
temperature dependence I/T2 is not clear yet. 

V. CONCLUSION 

Our ^^P NMR shift and susceptibihty data fitted reasonably well to the high temperature 
series expansion for 2D HAF model whereas the fitting to the ID model was not as good, 
especially in the T < 25 K regime. From the ^^P NMR shift analysis, the J/ks value 
was estimated to be about (18 ± 2) K. The large broadening of the ^^P NMR spectra at 
T < 5 K points towards a transition. However, no evidence of magnetic LRO was found 
in susceptibility and heat capacity measurements down to 2 K. Further experiments are 
required to really understand the detailed nature of this transition. ^■'^P NMR 1/Ti shows 
a good agreement with the theory of 2D HAF square lattice. The results reported in this 
paper thus suggest a variety of experiments and a need for a better theoretical understanding 
of quasi-low dimensional Heisenberg antiferromagnets. 




(8) 
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Figure Captions 

FIG. 1 (a) A schematic figure of the he plane in Na2CuP207 with [Cu(P04)2]^ linear 
chains propagating along c-direction indicated. A possible coupling path Cu-O-P-O-Cu is 
also indicated, (b) The arrangement of Cu and O in the 6c-plane is shown. Two planes 
formed by Cu^^ and Cu2^ ions are represented by thick and thin bonds respectively. 

FIG. 2 3ip NMR shift K vs temperature T for NaaCuPaOT. The solid line is fit to Eq. 
Hlin the temperature range, 15 K < T < 300 K, where Xspin is the high temperature series 
expansion for susceptibility of HAF square plane (Eq. Q). The dashed line is the simulated 
curve from Eq. [fusing the parameters obtained from the high-T fit. In the inset we have 
displayed the data at low temperatures on a logarithmic scale in order to show the deviation 
of experimental data from low-T series expansion. 

FIG. 3 3ip NMR shift K vs temperature T for NaaCuPaOr. The solid line is fit to Eq. 
E]in the temperature range, 5 K < T < 300 K taking Xspin for HAF chain (Eq. Ej). In the 
inset we have displayed the data on a logarithmic scale in order to show the deviation of 
experimental data from the theory around the broad maximum region. 
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FIG. 4 Low-T field sweep ^^P NMR spectra for Na2CuP207 are shown at different 
temperatures T around 5 K. It also shows the sudden change in line width and appearance 
of several distinct peaks. 

FIG. 5 Magnetic susceptibility {M/H) vs temperature T measured at H = 5 kG for 
NasCuPsOr. The solid line is best fit of the data to Eq. Elin the 15 K < T < 400 K range 
taking Xspin for HAF square plane (Eq. [T)). 

FIG. 6 Normalised specific heat Cp/{NAkB) of Na2CuP207 is displayed as a function 
of temperature T. The inset has d[Cp/ {NAkB)]/dT as a function of T showing an anomaly 
around 10 K. 

FIG. 7 The ^^P nuclear spin-lattice relaxation rate 1/Ti vs temperature T for 
Na2CuP207. The open circles are our experimental results and the solid line represents 
the simulated curve of Eq. [7| In the inset, the magnetization recoveries are plotted as a 
function of pulse separation t and the solid line is an exponential fit to Eq. |H1 

FIG. 8 The '^^P spin-spin relaxation rate I/T2 is plotted as a function of temperature T. 
In the inset, spin-echo decays are plotted as a function of t at three different temperatures 
for Na2CuP207. The solid lines show the fitting to a exponential function (Eq. |H1). 
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